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Dilute magnetic impurities in a disordered Fermi liquid are considered close to the Anderson metal-insulator 
transition (AMIT). Critical Power law coiTelations between electron wave functions at different energies in the 
vicinity of the AMIT result in the formation of pseudogaps of the local density of states . Magnetic impurities can 
remain unscreened at such sites. We determine the density of the resulting free magnetic moments in the zero 
temperature limit. While it is finite on the insulating side of the AMIT, it vanishes at the AMIT, and decays with 
a power law as function of the distance to the AMIT. Since the fluctuating spins of these free magnetic moments 
break the time reversal symmetry of the conduction electrons, we find a shift of the AMIT, and the appearance 
of a semimetal phase. The distribution function of the Kondo temperature Tk is derived at the AMIT, in the 
metallic phase and in the insulator phase. This allows us to find the quantum phase diagram in an external 
magnetic field B and at finite temperature T. We calculate the resulting magnetic susceptibility, the specific 
heat, and the spin relaxation rate as function of temperature. We find a phase diagram with finite temperature 
transitions between insulator, critical semimetal, and metal phases. These new types of phase transitions are 
caused by the interplay between Kondo screening and Anderson localization, with the latter being shifted by the 
appearance of the temperature-dependent spin-flip scattering rate. Accordingly, we name them Kondo-Anderson 
transitions (KATs). 



I. INTRODUCTION 



The non-Fermi liquid behavior of disordered electronic sys- 
tems such as the power-law divergence of the low-temperature 
magnetic susceptibility, can originate from a wide distribu- 
tion of the Kondo temperatures of magnetic impuri ties. Q^H The 
Kondo temperature Tk is exponentially dependent on the lo- 
cal exchange coupling J and thus on the hybridization of a 
magnetic impurity state with the conduction band. Since the 
hybridization is proportional to the overlap integral between a 
localized magnetic impurity orbital and the conduction band, 
it can be exponentially sensitive to microscopic variations of 
the position of a magnetic impurity. ' The exchange coupling 
also depends on all wave function amplitudes of the conduc- 
tion electrons at that position, and thus on the local density of 
states (LDOS) in the conduction band. The distribution of the 
LDOS depends strongly on the nonmagnetic disorder strength 
and is known to attain a wide log-normal distribution.'^l In 
the earliest approaches to this problem it was argued that the 
LDOS in the vicinity of the Fermi energy depends only slowly 
on energy. Therefore, the small Tk tail of the distribution 
should be directly connected to the distribution of the LDOS 
at the Fermi energy, which results in a wide distribution of Tk 
on the metallic side of the AMIT whose width increases as 
the AMIT is approached.^ However, since Tk is determined 
by an integral over all energies, the fluctuations of the LDOS 
are to some degree averaged out so that Tk may not vary as 
strongly. In the weak-disorder limit, explicit analytical cal- 
culations show that the width of the distribution of Kondo 
temperatures is finite due to correlations between wave func- 
tions at different energies. In disordered metals, these corre- 



lations are induced by the diffusion of conduction electrons. 
However, the resulting width was found to be small in three 
dimensions in the weak-disorder limit.l^ For strong disorder 
the distribution of Tk has been studied by means of the ex- 
act numerical diagonalization of finite systems to be wide and 
bimodaI.'S3its low-Tjf power- law tail has been argued to have 
a universal power. Similar bimodal distributions of Tk have 
been found in 2D disordered metals using the nonperturba- 
tive numerical renormalisation group and the quantum Monte 
Carlo methods.fifi' 

In this paper, we consider dilute magnetic impurities in a 
disordered Fermi liquid close to the Anderson metal-insulator 
transition (AMIT). We use information about the multifrac- 
tality of the critical wave functions at and in the vicinity of 
the AMIT,E2tS] in order to derive physical properties arising 
from the interaction between the conduction electron spins 
and the quantum spins of the magnetic impurities at the tran- 
sition point. In particular, we obtain analytical expressions for 
the distribution functions of the Kondo temperature Tk at and 
in the vicinity of the AMIT. 

At finite density, the magnetic moments become coupled 
by the indirect exchange coupling J^¥^ The latter can be 
calculated using the expression for the generalized RKKY 
coupling,'"* which is a function of the local exchange coupling 
J.GSlThe distribution of J^ and its competition with the Kondo 
screening will be studied in a subsequent article."' Here, we 
consider the dilute limit when the coupling between the mo- 
ments can still be disregarded. In the metallic phase the den- 
sity of free, unscreened magnetic moments npM is found to 
vanish exactly in the zero-temperature limit. In the insulating 
phase, npMiT = 0) is finite and increases with a power law 
as function of the distance to the AMIT. 



In Sec. In] we begin with a brief review of multifractal statis- 
tics. Wave functions in the vicinity of the AMIT are further- 
more power-law correlated in energy.EEHI As has been noted 
in Ref. 18 this has the surprising consequence that multifrac- 
tal eigenfunctions which are close in energy are likely to have 
their maximal intensity at the same positions in space. This 
positive correlation has been called stratification'^ and is op- 
posite to what is expected in the localized phase where states 
close in energy have their maximal intensity most likely at far 
apart positions in space. Another consequence is the open- 
ing of local pseudogaps at positions where the critical wave 
functions have vanishing intensity.^' At other sites the local 
intensity diverges with a power law. 



In Sec. Ill the Kondo problem in disordered systems is for- 
mulated. 
In Sec. 



IV the distribution of Kondo temperatures at the 



AMIT is derived and discussed. 

In Sec. [V]we extend the analysis to the metallic regime 
where the Fermi energy is located above the mobility edge. 
Although all states in the vicinity of the Fermi energy are 
extended, the multifractal nature of the eigenstates on length 
scales smaller than the correlation length ^ leads to fluctua- 
tions of the LDOS and modifies the magnetic properties of 
the system as we review in Sec. IV] 

In Sec. VI we extend the analysis to the insulating regime 
when the Fermi energy is located below the mobility edge and 
the wave functions are localized exponentially with a localiza- 
tion length ^c- Multifractal fluctuations still occur on length 
scales smaller than ^^ and the wave functions are power law 
correlated within a localization volume. It is therefore crucial 
to take these effects into account in order to get the correct 
distribution of Kondo temperatures in the insulating phase . 

In Sec. VII we summarise the results for the 2D system, 



where all states are localised. 



In Sec. VIII the quantum phase diagram of the Anderson 



metal-insulator transition in the presence of magnetic impu- 
rities is derived and plotted as function of the exchange cou- 
pling parameter J and disorder amplitude W. We establish 
the existence of a critical semimetal phase where both the 
correlation and the localization lengths are infinite within a 
finite interval of disorder amplitudes, and the conductivity is 
vanishing. 
In Sec. 



IX we consider how the quantum phase diagram 



changes in an external magnetic field, which couples via the 
Zeeman term to the magnetic impurities and via the orbital 
term to the conduction electrons. 

In Sec. |X]we present the results for the Non-Femi-liquid 
properties, in particular the magnetic susceptibility, the spe- 
cific heat, and the spin relaxation rate as functions of temper- 
ature, concentration of magnetic moments, and disorder am- 
plitude. 



In Sec. XI we study the consequences of the temperature 
dependence of the spin relaxation rate, which is caused by the 
Kondo effect. This leads to transitions between insulator, crit- 
ical semimetal, and metal phases at finite temperatures. These 
one may call, accordingly, Kondo-Anderson transitions. 



In Sec. Xn we provide our conclusions, comment on exper- 
imental realizations of these transitions, and discuss remain- 



ing open problems. 

In Appendix |A] we review the wave function correlations 
in the vicinity of the AMIT when one state is at the mobility 
edge. The joint distribution function of eigenfunction inten- 
sities is derived such that it matches the critical correlation 
functions. Next, the correlation function and the joint distri- 
bution functions are derived when both eigenstates are located 
away from the mobility edge of the AMIT. We also discuss 
and present results on higher moment correlation functions. 

In Appendix[B]we derive the function F{a, Tk) as defined 



inEq. (12 1 



II. 



MULTIFRACTALITY, LOCAL PSEUDOGAPS AND 
POWER LAW DIVERGENCIES 



The AMIT is known to be a 2nd order quantum phase tran- 
sition, where both the localization length and the correlation 
length diverge to infinity. At the AMIT the electrons are in 
a critical state, which is neither extended nor localized, but 
sparse, as seen in Fig. fTl where the intensity |?/'/(r)p at the 
AMIT is plotted. 



A. Multifractality 

This critical state can be characterized by the moments of 
eigenfunction intensities |?/;/(r)p, which scale as powers of 
the system linear size L, 
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where d is the spatial dimension. In a metal the powers Tq 
would be given by d{q — 1). Critical states are characterized 
by multifractal dimensions dq < d which may change with the 
power q of the moments. These are related to the exponents 
of the q-th moments hy Tq ^ dq{q — 1). The corresponding 
distribution function of the intensity is close to log-normal in 
good approximation,-^ 
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where a ~ — In |?/'i(r)p/lnL, r/ — 2{ao — d) and ao > d. 
The multifractal dimension dq is then related to ao by dq = 
d — q{aQ — d) for not too large q. At qc — ao/2/{aQ — d) 
there is a termination of Tq so that it remains constant Tq ~ Tq^ 
for q > qc\^ Throughout this article we assume the validity 
of this Gaussian distribution of a. In Fig. [Tithe local intensity 
is plotted for a critical state at the three-dimensional AMIT as 
obtained by exact diagonalisation of the Anderson tight bind- 
ing model with a box distribution of width W ~ 16.5t. Here, 
t is the hopping parameter. The energy is at £■ = 2t on a cubic 
lattice with spacing a, and linear lattice size size L = 100a. 
The coloring of the plotted intensity was done according to 
the amplitude of a = — In j-^p/ In L. Sites with higher inten- 
sity with a < 1.2 are so rare that their occurrence can not be 
resolved in this plot. All other sites whose intensities are not 
plotted correspond to lower intensity with a > 3.0. Thereby, 
about 80 % of the total state intensity is shown in Fig. fl] 
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critical energy Ek = Em is |7/'M(r)|^ = L~°'. This condi- 
tional intensity, relative to the intensity of that of an extended 
state, is obtained by averaging over the joint distribution func- 
tion Eq. M (for the derivation see Appendix A),*^ 



/„ - L''{\Mr)\^)lMr)\^=L- 
El — Em 



(5) 



where the power is given by 
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FIG. 1 : (Color online) The local intensity is plotted here for a critical 
state at the three-dimensional AMIT, as obtained by exact diagonal- 
isation of the Anderson tight binding model with W = 16. 5t, with 
hopping parameter t at energy _B = 2i on a cubic lattice with spacing 
a, and size L = 100a. The coloring of the plotted intensity was done 
according to red: a € [1.2, 1.8], green: a £ [1.8, 2.4], light blue: 



a € [2.4, 3.0], where is defined by a 



Inl 



/InL. Sites with 



higher intensity, a < 1.2, are so rare that their occurrence can not be 
resolved in this plot. All other sites whose intensities are not plotted 
correspond to lower intensity a > 3.0. Thereby, about 80 % of the 
total state intensity is shown. 



B. Critical Correlations 

The wave function intensity of a state at energy Ei at a 
certain coordinate r has a power-law correlation with the in- 
tensity at the mobihty edge energy Em with the power r]/d, 
which is related to a^ by 



77 = 2(ao - d), 



(3) 



(see Eq. ( |aT) in Appendix A).EG1I As noted in Ref. .22. this 
congelation is due to the fact that eigenstates in the vicinity of 
the AMIT have a similar multifractal envelope function (f>i{r) 
and can be written as tpi{r) = (j)i{r)xui{'r)- Here Xuii'r) de- 
notes a factor of the wave function which fluctuates indepen- 
dently from 0i(r) on short, microscopic length scales. 

Joint Distribution Function. One can take this correlation 
into account by considering the joint distribution function of 
two wave functions. In order to obtain the correlation function 
Eq. (JATJ and the distribution function of a single state, Eq. 
(|2]l, it should be of the form 

(4) 
where we obtained Kik ~ Max [\Ei — Ek\, A^^ ] /Ec, see Ap- 
pendix A for the derivation and definitions of aik- Here, ^; is 
the coiTelation length/ localization length of the state with en- 
ergy El on the metallic/insulating side of the transition, and 
A^ = D/^f. 

Conditional Intensity. Thus, we can derive the conditional 
intensity of a state at energy Ei given that the intensity at the 



This result is valid for \Ei — E]^j\ < Ec, where E^ is the en- 
ergy scale over which the critical correlations exist and giM = 
In [(E"/ — Em) / Ec\/ {d\n L). The average is done over the 
inte nsity [■(/;; (r)|^ using the conditional distribution function, 
Eq. ( |A7| . fixing the intensity L"" at the AMIT. TypicaUy, E^ 
is a fraction of the bandwidth L»: Ec ^ D/{2dln2d).^'' When 
El is located at a finite energy interval away from the mobility 
edge Em, the coefficient gij^i = In |(_E; — EM)/Ec\/{dlnL) 
vanishes for L — > 00. Close to Em the coefficient saturates: 
giAi\E,~>EM ^ ^1 ^iid Eq. (pjl reduces to L'^^", the local 
intensity at Em relative to the intensity of an extended state 
L^"^. Note that in Eq. (pb the average over the uncorrelated 
factor of the wave function Xui (r) has been done, which is 
{\Xui{v)?) = l. 

Local Pseudogaps. From Eq. (J5]l we see that at positions in 
space where the local wave function intensity at the mobility 
edge is small, corresponding to a larger than its typical value 
ao, the wave function intensity is suppressed within an energy 
range of order Ec around Em- Thereby local pseudogaps are 
formed with a power r^ = ""^"" with vanishing LDOS at the 
mobility edge, as shown in Fig. [2] 

Local Power Law Divergency. On the other hand, when the 
intensity at the mobility edge is larger than its typical value 
L^"", which corresponds to a < ckq, the LDOS is enhanced 
within an energy range of order Ec around Em, increasing as 
a power law when Ei approaches the mobility edge, as shown 
in Fig. |2] 



III. KONDO EFFECT IN A DISORDERED ELECTRON 
SYSTEM 

A. Kondo Impurity Hamiltonian 

Magnetic impurities are generally described by an Ander- 
son model where a localized level with energy e^ and on-site 
Coulomb repulsion U hybridizes with electrons in the con- 
duction band which is described by a Hamiltonian 7?o- This 
Hamiltonian may include the random potential of nonmag- 
netic impurities V{y). Using the eigenstates {V'n} and ener- 
gies {En} of iJg, with the corresponding one-particle density 
operator given by fin^a, the Anderson Hamiltonian is written 



\ 
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FIG. 2: (Color online) The conditional intensity /„ relative to the 
intensity of an extended state as function of the distance in energy to 
the mobility edge in units of the correlation energy Ec- The exponent 
a, which is related to the intensity L^" at Em, takes the values 
a — 2,3,3.5,4,6 in the sequence of decreasing dashing. Local 
pseudo gaps are seen for a > ao, while local power law divergences 
occur for a < qq- Here, we set d = 3 and ao = 4. 



B. Kondo Temperature 

As we have shown previously with the numerical renormal- 
ization group and the quantum Monte Carlo methods,'" the 
distribution of Kondo temperatures is in very good qualitative 
agreement with the one obtained from the one-loop equation 
of Nagaoka and Suhl.^'' Only, its average needs to be rescaled, 
in order to account for a shift of the distribution to ward s larger 
values of Tr- due to the higher-loop corrections 112122] There- 
fore, we will use here the one-loop equation to calculate Tk, 
as given by 
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F{{ai},TK). 



(10) 
where j — J /D, A = D /N is the mean level spacing and 
ei — El ~ En. N is the number of states in the sample of 
volume V = L''-. This defines the Kondo-temperature in terms 
of the local intensities at all energies e; in the sample. 



C. Distribution Function of the Kondo temperature 



Ha = '^Enfin.a + ed'^hd^a + Und,+nd- 

n^a (7 

+ ^{tndC+„d„^Y{.c), (7) 



where fid.cr is the density operator of the impurity level. 
The hybridization amplitude tnd is proportional to the eigen- 
function amplitude at the position of the magnetic impu- 
rity V'n(O) ^nd to the localized orbital amplitude (t>d{^)'- 
tnd = iV',*(0)<^d(0). One can employ the Schrieffer- Wolff 
transformation,-^ ^^ formulated in terms of eigenstates ijjn to 
take into account double occupancy up to second order in 
^nd- The result is an s-d contact Hamiltonian with exchange 
couplings given by 



f* i , 
''nd ^ri'a 



1 



1 



U + ed — En' En 



fid 



(8) 



and an additional potential scattering term with amplitude 



''nd'^n'd 



1 
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U + ed — En' En — ed 



(9) 



Note that in the symmetric Anderson model, Knn' vanishes 
for all n,n' when \En.n' — Ep\ ^ U, since, in this case 
td = Ep — U/2, where Ep is the Fermi energy. For arbitrary 
ed wave functions with small amplitude at the position of the 
magnetic impurity are hardly modified by this potential scat- 
tering term since Knn' ^ V'i*(0)V'n'(0). Hence, we will first 
retain only the exchange couplings J„„' ~ Ji/),* (0)?A„'(0), 
with J ^ t'^/U, and leave the discussion of possible effects 
of finite Kn n' to Sec. IXII 



Thus, one can derive the distribution function of 
Tk, when the distribution function of all intensities 

P{ai, ..., ai, ...,aAr) is known, by solving 

P{TK)^Jl[daiP{{ai})Sil-Fi{ai}^TK))\^l 

(11) 
where F{{ai},TK) is defined by Eq. (fTol. Note that 



F{{ai},TK) is always a decreasing function of Tk, so that 

I dF I _ dF 
I dTK I ~ dTK- 



IV. KONDO EFFECT AT THE 
ANDERSON-METAL-INSULATOR TRANSITION 

Conditional Average. In a first attempt to get its distri- 
bution function we calculate the Kondo temperature Tk for a 
given intensity i^~" at the Fermi energy, and integrate over 
all other intensities with the conditional distribution func- 
tion Eq. (A7i for fixed a. Thereby we find TK{a) as a 



function of a, using that the conditional intensity of state I, 

(l'/'/(r)H|^„(r)p= 



IS given 



by Eq. (Isll, 



1=J 



A 



E 



tanh 



(I 
2Tk 



= F[a,TK], (12) 



where the summation over I is restricted to energies within 
the energy interval of the correlation Energy Ec around the 



mobility edge 



l28l 



We note that Eq. ( 12 1 defines the Kondo temperature in a 



system with pseudogaps of power r^, Eq. (J6]l, in the local den- 
sity of states shown in Fig. 2, when this power is positive.l^Sj 
Therefore, the Kondo temperature is reduced at such sites, 
since the magnetic impurities are placed in the locally sup- 
pressed LDOS of the conduction electrons. 



On the other hand, at sites where the power r^ is negative 
the opposite occurs and Tk is enhanced. The Kondo prob- 
lem with a power-law divergence in the density of states was 
studied in Ref. 1^^ finding an enhancement of Tk towards the 
strong coupHng fixed point given by J. Since a is distributed 
over all values from to infinity according to Eq. Q, we find 
by solving Eq. (12 1 that Tx{a) is distributed, accordingly. 



We note that Eq. (10 1 is averaged over the uncorrected part 
of the wave function, x„i(r). Since these short scale fluctua- 
tions are uncorrelated in energy, they are averaged out by the 
summation over energy levels I in Eq. ([12]), and gives rise 
to small fluctuations of Tk of order of 1/N, only. Therefore 
their effect on the distribution of Tk is negligible in the ther- 
modynamic limit. 

We can now proceed and solve Eq. ( 12 1 analytically in 
various limits. The Kondo temperature for a given intensity 
|'(/'„(r)p — L^" at the Fermi energy is for a > ao — d found 
to be given by 



a — ckq 1 

2a. — 7] 



Eq. 



r.Tn 



(13) 

takes into 
deriving it 



Tk 
Eo ~ 

Here, j = J/D and .„ - ^=^2+3' 
account the critical correlations at the AMIT." 
(see Appendix IB] for details), we approximate tanh x ^ x 
for X < 1, and tanh a; « 1 for a; > 1, yielding two terms. 
One includes the integral over all energies within a window 
of width 2Tk around the Fermi energy. The other one is over 
all larger energies in the conduction bands. In order that the 
resulting equation is valid for all a it is important to keep both 
terms. 

For the typical value a = ao we recover Tk of a clean 
system in the one-loop approximation, namely. 



TK{a = ao) ^ £^cexp(-l/j) -- T\ 



(0) 
K ■ 



(14) 



From Eq. ( 13 i we see that at particular sites where the wave 
function amplitude is large, corresponding to a < a^, Tk 
becomes enhanced. In particular, for a wave function intensity 
comparable to the one of a metallic state, a = d, we find that 



TK{a = d) ^ j 
This is larger than T^ ' when j < 1 



2/v/d 



(15) 



Out of the N = L'^ atomic sites in the system, magnetic mo- 
ments remain free at all temperatures if placed on one of NpM 
sites, where a sufficiently strong local pseudo gap is devel- 
oped. Thus, for small J there can be a macroscopic number 
of such sites, although their density npM — N^m/N, 



UFuiT ^ OK) ^ L-'' ^nu, 



(17) 



vanishes, npM — > for L 



B. Distribution Function of Kondo Temperature. 

1. Tk ^0 Limit of P{Tk) 

In deriving the conditional intensity, Eq. (J5]l, and inserting 
it in Eq. ([12]) we took into account the correlations of all states 
to the intensity at the Fermi energy as characterised by local 
pseudogaps and power law divergencies. Since their power Tq, 
is distributed, we can now obtain the distribution of the Kondo 



temperature Tk by solving Eq. ( 13 i for a(T/f ) and inserting 
it in P{a). To this end, we use the Fourier representation 
of the delta-function. Then, we can expand in F[{a;}, Tk] 
and perform the average, keeping ai = a at the Fermi en- 
ergy Ep ~ Em fixed. Next, we use the conditional pair ap- 
proximation introduced above. If we impose the condition 
F[{q;;}, Tk] = 1 to obtain Tk for given a and insert this in 
P{a{TK)) this yields in the limit of small Tk -^ A, 



P°{Tk ^ 0) 
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(18) 



Forj = 0.25 the power of the Tk -^ Otailis/3 = 1-j = 0.75 
in exact agreement with the numerical result in 3 dimensions 
reported in Ref. 8 However, we find that its weight is vanish- 
ing with a power of the system size L. Note that the number 
of sites used in Ref. 8 is N = 2197, yielding a level spacing 
A/T^ ~ 0.025, so that Eq. ([is]) indeed can explain the tail 
of the distribution Tk < A displayed in their Fig. 3. Clearly, 
at larger Tk the fluctuations of the wave function intensities 
at energies away from the AMIT are important. These can 
strongly change Tk and the width of its distribution as we 
find in the next subsection. 



A. Density of Free Magnetic Moments at Zero Temperature. 

At sites where the wave function amplitude is small, corre- 
sponding to a > ao, Tk is suppressed due to the appearance 
of local pseudogaps. The Kondo temperature in the presence 
of pseudogaps of power r is well known to vanish when the 
exchange coupling does not exceed a critical value jc = r^ 
Since the power of the local pseudogaps depend on a, Eq. ([6]), 
that critical value jc{ct) depends on a as well. Accordingly, 
magnetic moments remain unscreened even at T = OK, when 
a exceeds the critical value 



apM = ao + dj. 



(16) 



2. P{TK)forTK> A 

In order to proceed in the calculation of P{Tk), at Tk ex- 
ceeding the level spacing A we need to take into account the 
fluctuations of the intensities at all energies. To this end, let 
us first rewrite Eq. ([TT) as 
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(19) 



where (...)q, denotes the average over all intensities, except 
the one at the Fermi energy which is fixed to L^". Next, we 



use the Fourier representation of the delta-function to get 



P{Tk) = 



daP{ 



dt- 



x{exp{-itF[TK]))a- (20) 



Now we can expand F[{ai},TK], around F[a,TK] = {F)a 
to take into account the fluctuations of all intensities. Expand- 
ing in 6F = F[{ai},TK] — F[a,TK], and performing the 
integral over t, we find 



PiTK) 



da , ,dF\a,TK\ (i--F[°.^gi)^ 



where F\a, Tk] = {F)a, is given by Eq. (B2 1 and F is de- 
fined by 



F^ = (F2)„ - {F)l, 



(22) 



where (...)„ denotes the average over all ai, keeping only a 
at the Fermi energy fixed. We find, that for A < Tk ^ Ec, 
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(23) 



where we determined numerically, ci « 7.51 and C2 ~ 9.60. 
Note that F vanishes in the hmit Tk -^ A, T{Tk ^ A) = 0, 
since then only energy levels in a range of order A around 
the Fermi energy contribute, whose correlations are taken into 
account correctly by Fa, already. Thus in this limit, the con- 
dition 1 = F{a, Tk) is imposed exactly, and we recover the 
tail of the distribution, Eq. (18 1, diverging with the power 

At larger Tk > A, F has a finite value and decays for A < 
Tk <S Ec with the power rj/d. P{a) is peaked ata = oq with 
a width that scales with the system size as l/y/hiL. Thus, for 
L — > oo, a = ao is imposed in Eq. (21 1 for any finite Tk- 
Thus we can substitute 
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(0) 



and find that the distribution diverges at A < Tk ^ T\^ as 



P{Tk) 




(25) 



with the power (3 — 1 — ri/d/2. In d = 3 dimensions, with 
ao = 4, the power is I3d=j, = 2/3 which is smaller than the 
one obtained numerically in Ref. |8] /3 ~ 0.75. We note, 
that there is a noticeable deviation towards smaller powers for 
Tk > A « 0.03 in the Fig. 3 of Ref. [3 



a. 




FIG. 3: The distribution of a defined in the metallic regime by Eq. 
p.6\ , as obtained by exact numerical diagonalisation of a 3D sample 
of size L"^ = 128'' in units of the grid cell volume a^. The energy 
of that state is approximately E = 0, and a box distribution of un- 
correlated disorder potential with W — 15t is taken. The analytical 
expression Eq. \27\ is plotted for comparison as the solid line for a 
correlation length ^ = 22a (taken from the numerical results of Ref. 

unii. 



V. KONDO EFFECT IN THE METAL PHASE 

A. Multifractality in the Metallic Phase 

In the metallic regime all wave functions are extended 
and their intensities scale with the inverse system volume, 
I^P ~ L""*. On length scales smaller than the correla- 
tion length ^, multifractal fluctuations of the wave function 
intensity occur as lo ng as ^ is larger than the microscopic 
length scale Ac.'^^" As pointed out in Refs. I20I22I the mo- 
ments of the intensity do scale with the correlation length ^ as 
L'^i{\^\^i) - ^(<i-d,)(g-i). Therefore, in the metallic regime 
we define a as 



L''\^i{vr=c^ 



d—ai 



(26) 



where f; is the correlation length of state /. Notice that this 
definition of a crosses over to the one we used above in the 
critical regime, where ^ diverges and is replaced by the system 
size L, when L < ^. It has to a good approximation still the 
Gaussian distribution. 



P{ai) r^ exp[-ln^i 



2,7 



(27) 



where its width scales with the logarithm of the correlation 
length ^i. This is confirmed in Fig. [3] where we plot Eq. 
( 27 1 together with the numerical result, as obtained from exact 
diagonalisation. 
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FIG. 4: (Color online) The conditional intensity la at energy Ei, 
given that a state at the energy Ek has intensity L''\ilik{f)\^ = Cf""' 
relative to the intensity of an extended state. It is plotted as func- 
tion of the distance of energy Ei from the mobility edge energy Em, 
in units of the correlation energy Ec- The exponent a at the en- 
ergy Ek, takes the value a = 4.5. The energy Ek is varied from 
the mobility edge into the metallic regime with {Ek — Em)/ Ec = 
0, 0.2, 0.3, 0.4, 0.7, from left to right. Instead of a local pseudo gap, 
one sees an increasingly shallow depression. Here we set d = 3 and 
Qfo = 4. Note that it varies slowly in a range Aj^ around the Fermi 
energy. 
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FIG. 5: (Color online) The distribution of Kondo temperatures Tk 



r.{0) 



in units of T]^ in the metallic phase, Eq. 1 33 i, is plotted as function 



of the distance to the mobility edge, Ep — Em in units of Ec for an 
exchange coupling j = 1/5. 



B. Intensity Correlations in the Metal 



There are still power-law correlations in energy between 



wave function intensities as given by Eq. (A3 1. Averaging 
with the conditional distribution function Eq. ( |A4[ ) in the 
metallic regime, we find the conditional intensity in the metal- 
lic regime: given that a state at the energy E^ has intensity 
L''|-0fe(r)p — Cj. ", a state at energy Ei at position r has on 
average the intensity 



Uiu^k) = (il^/WHo = kI '"" ^ "''^" , (28) 



where Kik = Max {\Ei- Ek \ , Min(A^, , A^ J} /E, and 

Aj, = E.iaJ^i)", (29) 

is the mean level spacing of a sample of finite size S^i (for the 
derivation, see Appendix [A|. Thus, the intensity at the energy 
El has still a dip when a > ao, although the LDOS at energy 
El —> Ek is no longer suppressed to zero but rather to a finite 
value (see Fig. l4]i, given by, 



C. Distribution of Kondo Temperatures 

Since the intensity is finite at all sites in the metallic regime 
we expect that the magnetic impurity spin is at low temper- 
atures always screened. Therefore, the Kondo temperature 
does not vanish in the metallic regime at any site. According 
to Eq. ( [30| the local intensity can be substantially suppressed 
or enhanced, depending on the random value of a. Its distri- 
bution, as given by Eq. (27 1, has a finite width. On the other 
hand, the second moment of F does saturate to a finite value 
whenTft: < A^: 



r^lT.<A,^Ai(-^«^ 



Er, 



-v/d 



(31) 



where ci k, 7.51. Inserting both Eqs. ( 31|27 i into the ex- 
pression for P{Tk), Eq. (21 1 we find that the probability 
to find a small Kondo temperature such that Tr- < Aj = 

\{Ef — Em)/Ec\'"'', is decaying to zero as 



P(Tj 



K 



In-^ 



Tk 



.(32) 



{L'mr)]'). 



Ei—^Ek 



i 



(d-a) 



Ek - E, 



M 



Er^ 



iy{a — d) 



(30) 

which is slowly varying with energy in an interval of order 
Aj^ around the Fermi energy. 



Thus, in the metallic phase Ep > Em there are no free mag- 
netic moments and the low T/v'-tail terminates at Tk ~ A^ 
as seen in Fig. (jSl). For Tk > A^ a power law tail of the 
distribution can still be observed, as seen in Fig. (Bl. There, 
we plot P{Tk) as obtained by Eq. (25 i and substituting Tk 



by Max(TA', A^) in the expression for F, Eq. (23 1, 



PiTx) 



Er. 



Tk 



exp ■ 



J_ r Max{TK,A^) 
2ci I E, 
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VI. KONDO EFFECT IN THE INSULATOR 

In the insulating regime each locaHzed state is restricted to 
a volume of the order of S,f = £,^{Ei), the localization vol- 
ume. On length scales smaller than the localization length ^; 
there are still multifractal fluctuations. In addition, the wave 
functions in the insulating regime are power-law correlated 
close to the AMIT.'^ Neglecting a small logarithmic enhance- 
ment which occurs at energy spacings smaller than the local 
level spacing A^ — S^^'^D^ we can get the distribution of 
Tk simply by using the results obtained at the critical point 
and replacing the system size L by the localization length ^c- 
Defining 



\M^)\^ = ^i 



(34) 



we thus find that the distribution function of a within a lo- 
calization volume is the same as if we had considered the 
distribution at the AMIT in a finite volume of order i^;*^. We 
note that the probability to find a state at energy Ei inside 
the localization volume is decaying with the system volume 
L'^ as ^f /L'^. Outside the localization volume the intensity 
decays exponentially, |-0;(r)P ~ exp[— 2r/^] corresponding 
to a = 2r/(^;ln^;). Since most sites in the sample are a 
distance r ^ L — > oo away from the localization volume one 
finds a ^ L — > oo at most sites of the sample. The correlation 
function between two wave functions at different energies Ei 



and Ek is still given by Eq. ( Al i where A^, is now the local 



level spacing at energy Ei. Accordingly, the joint distribution 
function has the form given in Eq. (A4i. However, the dif- 



ference to the metallic regime is that there are only discrete 
energy levels which are separated by the local level spacing 
A J, . This difference is important, especially when calculating 
the Kondo temperature, since the hard gap A^, in the insulat- 
ing regime cuts off the Kondo renormalization flow at small 
energies.-^' 

Thus, we can conclude that there is a finite density of free 
magnetic moments tifm, which remain unscreened. There- 
fore, we need to subtract this density from the distribution 



function of T^ given by Eq. ( 2 1 1. The width F is given by Eq, 



(23 I when substituting A by the level spacing in a localization 



volume, A^ for Tk > Aj, 

r^ ^ f ^^%,^'^\ .{^-^^)-^'\ (35) 

where Ci w 7.51. For Tk < A^, we find that it vanishes, F — 
0, and the condition 1 — F^ is enforced exactly. Therefore, we 



get from Eq. (21 1 the low Tk tail in the localised regime as 



P{Tk) 



npAiiT ^ 0) 



flAI 



E^ 
Tk 



i-j 



e 



-^(djf 



(36) 

Thus, we find that the distribution diverges with the power 
/3 — 1—j. This is in full agreement with the numerical results 
of Ref. 8 , where a power 0.75 has been obtained for j = 0.25 
for a wide range of disorder amplitude W. Also, the increase 
of the weight of the power law tail with disorder strength W 
is in good agreement with the numerical results. 

The finite density of free magnetic moments is found to 
be,'^ 

7ipm(t = ox) = n^r^^*^' 



um 



V Wc 



(37) 



which decays to zero as a power law when the disorder ampli- 
tude W approaches the AMIT at Wc- It converges to the total 
density of magnetic moments um far away from the mobility 
edge, Em — Ep — > Ec- Thus, according to this expression 
all magnetic moments should be free in the strongly localized 
regime where ^ — > Oc. 



VII. KONDO EFFECT IN 2D ANDERSON INSULATORS 

We can apply this analysis also to two-dimensional disor- 
dered electron systems, where all states are localized. The 
Kondo effect in such systems has been studied numerically 
based on the 1-loop equation,^ and with nonperturbative 
methods in Ref. [TO] where both the distribution of Kondo tem- 
peratures and the density of free magnetic moments have been 
obtained. The two-dimensional localization length in the ab- 
sence of a magnetic field is known to depend exponentially on 
the disorder strength, £^2D — 9 exp(7rg), where g = EpT. The 
scattering rate 1/r is related to the disorder amplitude W as 
1/r = ttW^/6D. There are weak wave function correlations 
in two dimensions which are logarithmic, with an amplitude 
of order 1/g. For weak disorder, 5 ^ 1 we can rewrite this 
correlation as an effective power law with power 



772D = 2/Trg. 



(38) 



The correlation energy in 2D is of the order of the elastic scat- 
tering rate, Ec 20 ^ 1/t. Thus, for systems whose size L 
is smaller than the localization length ^2D, the 2D system be- 
haves like a critical system with ao defined by 772D = 2/7rg = 
Qo — 2, or 



ao 



2 

ng' 



(39) 



r(i) 



There is a critical exchange coupling Jc above which there 
is no more than one free magnetic moment in the whole 
sample."^ Substituting 7720, we find, that 



j^y 



D 

3E^ 



W. 



(40) 
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This is in good agreement with Ref. [TO] where Jc has 
been determined numerically for a 2D disordered system, and 
found to increase linearly with disorder amplitude W. There, 
and in Ref. 8 the density of free moments have been obtained 
which we can now compare with our analytical expression. 



npuiT = OK) = nM^2D^" = "M(.gcxp(7r.g))" 



(41) 



with .9 = 6/{ttW^)Ep/D. By substituting £,2D intoEq. (|36i 
and inserting the parameters used in Ref. 10 we compared the 
analytical distribution of Tk with the numerical results and 
found a good qualitative agreement (not shown). We note that 
there is a low Tk power-law tail divergence where the power 
is given hy /3 = 1 — j. 



Vin. ZERO TEMPERATURE QUANTUM PHASE 
DIAGRAM 



It is well known that the scattering of conduction electrons 
by magnetic impurities can lead to the relaxation of the con- 
duction electron spin, and thereby the loss of the electron 
phase coherence. At finite temperature this leads to a suppres- 
sion of the quantum corrections to the conductance, the so- 
called weak localization corrections. In the low-temperature 
limit phase coherence is restored, but the magnetic scattering 
may still break the time-reversal symmetry of the conduction 
electrons similarly to an external magnetic field. The break- 
ing of time-reversal symmetry is known to weaken Anderson 
localization and thereby the localization length becomes en- 
hanced. In systems with an Anderson metal-insulator transi- 
tion the transition is shifted toward str onger disorder ampli- 
tudes W and lower electron density n. Thus, the sym- 

metry class changes, shifting the AMIT from the orthogonal 
symmetry class (time-reversal symmetric) to the unitary sym- 
metry class (broken time-reversal symmetry).'32llSl 

In the presence of an external magnetic field this change of 
the symmetry class of the conduction electrons from orthog- 
onal to unitary is governed by the parameter Xb = £,^/Ib' 
where Ib is the magnetic length. Therefore, the spin scatter- 
ing rate due to magnetic impurities 1/ts is expected to enter 
through the symmetry parameter Xg = S,^ /D^Tg, where D^ 
is the diffusion constant and ^ is the correlation (localization) 
length on the metallic (insulating) side of the AMIT. ^^ When 
Xg >1, the electron spin relaxes before it can cover the area 
limited by ^ and the system is in the unitary regime. One can 
then study the crossover of the mobility edge through a scal- 
ing Ansatz for the conductivity on the metallic side, as done 
in Ref. [32] in the case of a magnetic field. Following this 
approach, using the spin scattering rate l/r^, we get 



'^(IM) - ^/(Xs 



(42) 



The conductivity then goes to zero at the critical disorder 

W^'{1/ts) as CT(l/r,) ~ {W^'il/Ts) - Wf , where the in- 
dex M indicates that this is the critical disorder strength when 
it is approached from the metallic side, W < W^'^ (I/ts). 



Coming from insulating side, we need instea d to a pply the 
scaling Ansatz to the dielectric susceptibility t^^ P^P?! 



Xe{l/rs)=e9iXs). 



(43) 



This diverges at the critical disorder W^{1/ts) as x(1/ts) ^ 
{W — W^{1/ts))^^'^, where the index / indicates that this 
is the critical disorder strength when it is approached from 
the insulating side, W > WI{1/ts). While in a magjietic 
field these two critical points are found to coincide, ^-^ we 
will see below that, in the presence of magnetic impurities, 
W^'''{1/ts) and T4^/(l/r^,) can be different due to the Kondo 
effect. 

When a finite concentration of classical magnetic impuri- 
ties um with spin S is present, the magnetic relaxation rate at 
zero temperature is given by i/^-ciassicai _ 2TmMS'^pp{eF), 
where p{tF) is the density of states at the Fermi energy. How- 
ever, the quantum mechanical nature of the impurity spins af- 
fects this rate in several ways: First, its magnitude is enhanced 
since the quantum mechanical eigenvalue of the square of the 
spin is S{S + 1). This results for S = 1/2 in a factor of 
3 enhancement. Secondly, the Kondo effect tends to screen 
the impurity spin leading to a vanishing spin relaxation rate at 
zero temperature when magnetic impurities are dilute. How- 
ever, at finite temperature the Kondo correlation can instead 
enhance the spin relaxation rate with a maximum at Tk ■ This 
effect has been observed in weak-localization experiments as 
a pl ateau in the temperature dependence of the dephasing 
time!^^ Recently , its full temperature dependence was ob- 
tained numerically.EE2l A good agreement with the numerical 
results can be obtained through the approximate expression,!^ 



.(0) 



(T) 



■Kn,nS{S +1) J 2 



-tt'^S{S + 1) 



Tk 
T 



T 
Tk 

2 



(,44) 



with /? = 0.2, as obtained numerically.™ The temperature de- 
pendence scales with Tk- Note that in the low-temperature 
limit the spin relaxation rate vanishes as T'^/T^, similarly to 
the inelastic scattering rate in a Fermi liquid and in agree- 
ment with Nozieres' renormalized Fermi liquid theory of di- 
lute Kondo systems.Eil 



A. Metal Phase 

In the zero-temperature limit, coming from the metallic 
side of the AMIT, the Kondo screening results in the vanish- 
ing of the spin relaxation rate, 1/ts{T = 0) = 0.-- There- 
fore, Xg = and the AMIT occurs at the orthogonal criti- 
cal value for time-reversal symmetric systems, W^. For the 
Anderson tight-binding model in three dimensions, one finds 
Wp /t = 16.5 ± 0.02, with the orthogonal critical exponent 
i^o = 1.57 ±0.02 and the correlation exponent rj/do — 
0.56 ± 0.02.'^^"*. Most recently with a more accurate mul- 
tifractal scaling method W° /t = 16.530(16.524, 16.536), 
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vo = 1.590(1.579,1.602), ao = 4.048(4.045,4.050) and 
rio = 1.763(1.792, 1.727) was obtainecP. 

Looking now at the transition point we can conclude from 
the analytical results of Sec. VI namely Eq. ( 37 1, that at 



the AMIT the density of free magnetic moments vanishes. 
However, there can remain a macroscopic number of free mo- 
ments. We find that from the Nm magnetic moments in the 
sample, at least \/Nm of them remain free when the exchange 
coupling does not exceed the value jc = yil/d. 



Log PpM 




FIG. 6: (Color online) The fraction of free magnetic moments Pfm 
at T = QK, Eq. \?>1) , in a three-dimensional disordered metal as 
function of the exchange coupling J (in units of the band width D) 
and disorder strength W (in units of the critical value W^). Critical 
correlations result in a finite Pfm even for large J > J^ (dashed 
line]). For J < J* , Eq. \5Q\ there is a critical region for disorder 
amplitudes W^ < W < Wc{J), where Wc{j = J/D) is given by 
Eq. J49}. This figure was previously published by the authors in Ref. 



In Section IVllit was found that on the insulating side of the 
AMIT the LDOS at the position of unscreened magnetic mo- 
ments scales as p{eF, r) ^ ^d-a{j)^ where apj^i = aQ + dj is 
the critical value. Therefore, for a > apM the magnetic mo- 



ments remain free. According to Eq. (45 1 the spin relaxation 
rate due to the free moments depends itself on the localization 
length as 



1 



rFM 



(T = 0) 






p{<^f) 



^ 



(46) 



where the density of free moments is given by Eq. (37i 
and, depends on the localization length ^ as ripM = 

d -2 

'T'Af (^/oc)^^^^"' ■ Following Ref. 32, we get the critical 
disorder amplitude when the argument of the scaling func- 
tion g in Eq. ( [43] l is of order unity, yielding the condition 
Dts{J) = £.[Wj{J) - W^]'^ as function of the exchange 
coupling J. This condition is valid as long as the deviation 
from W^ is small. For larger deviations it will converge to 
its unitary value W^ . The exact analytical form cannot be 
obtained from this phenomenological scaling approach. We 
note that, in contrast to the case of an external magnetic field, 
as considered in Ref 32 according to Eq. ( [47] ), Tg depends 
on the localization length ^ as 



1 / fX-winU+n/df 
(T = 0) = ^^ ' ^''^ 



rFM 



where 



-^ = 2TTnMS{S + l)fp{tF) 



(48) 



B. Insulator Phase 

On the insulating side of the AMIT the density of free mag- 
netic moments is finite as given by Eq. (37 1. Therefore, the 



time reversal symmetry and the spin symmetry of the conduc- 
tion electrons is broken in proportion to the spin relaxation 
rate l/r^. This leads to a shift of the AMIT from VF° to 
WI{1/ts), which we will determine in the following. 

Before proceeding we note that in a disordered system the 
spin relaxation rate is proportional to the local density of states 
p{E^ r) at position r, since it is given by 



1 

Ts{E,T,y) 



^rrr ^ ^ 



p{E,v) 



PiE,r) 



'tllmGd + TipttGl 



r^\E,Ty 



(45) 



with the propagator of a localized d-level in the Anderson 
model Gd-"^^ Here, I/t] ' is the spin scattering rate in a clean 



system as given by Eq. (44 1. Thus, Eq. (45 1 leads us to con 



elude that the spin scattering rate not only depends on the 
ratio T/Tk but also explicitly depends on the LDOS. Since 
both Tk and p{E, r) are randomly distributed, \/ts{E, T, r) 
is distributed as well. 



Thus, we finally get the shift of the critical disorder as function 
of J as 



W^{j)^W^ + W', 



D, 



«(j) 



(49) 



where l/n{j) — v\2. — |-(j + rj/dY]. This result is vahd 

for small deviations from W^ . For larger deviations it will 
approach the unitary value W^ in a still unknown form. We 
see that for j = j* where 



J 



2VV V 
d d' 



(50) 



which gives in li = 3, j* w 0.276, the exponent k{j) di- 
verges and W^ approaches its orthogonal value W^, see Fig. 
p] Also, for larger values of J, it will stay at W^ as a conse- 
quence of the increase of Kondo screening with the exchange 
coupling J. 



C. Critical Semlmetal Phase. 

For smaller exchange couplings j < j* a paradoxical situa- 
tion appears: The position of the critical point W^ depends on 
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the direction from which the AMIT is approached. Therefore, 
for intermediate disorder strengths, W^ < W < W^{j) there 
exists a critical region. Accordingly, the mobility edge is ex- 
tended to a critical band whose width is a function of j. The 
resulting zero-temperature quantum phase diagram is shown 
in Fig. 1^ According to Eq. ( [49| with Eq. BSj, the width of 



the semimetal phase W^ (J) ~ W^c 



'Af 



power of the density of magnetic impurities um- 



IX. KONDO-ANDERSON TRANSITIONS IN A MAGNETIC 

FIELD 

A Zeeman magnetic field polarizes the free magnetic mo- 
ments. Thereby, their contribution to the spin re laxation rate 
becomes diminished by the magnetic field.EEiMZlon the other 
hand, the Kondo singlet which Kondo screened magnetic mo- 
ments form with the conduction electrons is partially broken 
up by the Zeeman field. Thus, these magnetic moments con- 
tribute a spin relaxation rate which is increasing with the Zee- 
man field. Finally, an orbital magnetic field breaks the time 
reversal symmetry and therefore also results in a shift of the 
AMIT towards the unitary limit. It is therefore an intriguing 
problem how these competing magnetic field effects combine 
to change the quantum phase diagram.!^ 

The magnetic field dependence of the spin relaxation rate 
from magnetic impurities in a metal at finite temperature was 
calculated in Ref. 147] The magnetic field polarises the mag- 
netic impurity spins due to the Zeeman interaction. 



H, 



-7. 



:B y ^ Si, 



(51) 



where 7s = .9s Mb is the gyromagnetic ratio of the mag- 
netic impurities, gs their g-factor and fiB is Bohr's magne- 
ton. Here, the magnetic field B is taken to point in the z- 
direction. The spin relaxation rate from free magnetic im- 
purities is found to be exponentially suppressed according to 
1/ts ^ ex-p{~'^sB\Sz\/kBT). Thus, in the zero temperature 
limit all free moments are expected to become polarized by an 
arbitrarily small magnetic field, and their contribution to the 
spin relaxation is vanishing identically. 

Magnetic impurities which are screened with a finite Kondo 
temperature Tk, however, have without magnetic field a spin 
relaxation rate Eq. ( [44| which vanishes in the low temper- 
ature limit. Applying a magnetic field the Kondo singlet is 
partially broken up and a fi nite spin relaxation rate appears, 
which scales with Tk ag^ED 

riM hsB\S,\)^ 



-(Tk) 



np 



rp2 



(52) 



To get the total spin relaxation rate, we integrate over the dis- 
tribution of Kondo temperatures P{Tx)- We note that the 
contribution from magnetic impurities with a small Kondo 
temperature, Tk < gslJ,BB\Sz\ vanishes since these spins 
become polarised. The ones with larger Kondo temperatures 
then yield the spin relaxation rate, 

dTKP{TK)^^ ^^2 • (53) 



1 



For small magnetic fields, ^sB\Sz\ ^ Tk\ the main contri- 
bution comes from the low T/^-tail of the distribution. 



p™H^)"r^^^ 



increases with a Thus, we get 



j_ ^ riM dj |^(j+^/rf)2 / jsB\Sz\ '^^ 

Ts 7Tp2-j \ Ec 



(54) 



(55) 



Setting Xs — ^ jDf^Ts = 1, we find that the Zeeman field 
shifts the critical disorder to 



W,{B) = Wlij) + W°CM 



lsB\Sz 



j«(i) 



(56) 



where 1/k(j) = 2^.(1 - "^{j + v/df) and cm = 

(djriM /{{'2, — j)npDe)) . Thus, the transition between 
critical semimetal and insulator is shifted in a magnetic field 
according to Eq. (|56]l as plotted in Fig. (|7]). 

The orbital magnetic field is known"*^ to shift Wc to 



WciB) = Wy + W^{-KeB/h) 



Nl/(2<.) 



(57) 



This determines the transition line between metal and 
semimetal, since, the Zeeman field contributes a slower de- 
pendence on B, coming from the metal side of the transition 

we find W^{B) = W^ + W^ i^^^^f^) ^'^ ■ 

For the transition between semimetal and insulator, we can 
conclude that for 

J < Jz = v/d {2y/d + 1 + dyr^/d - 1 - 2/d) , (58) 

the shift of Wc is dominated by the Zeeman field over the 
orbital magnetic field. For, d — 3, rj/d — 2/3 one finds, 
jz = 0.185, so that the Zeeman field effect dominates for 
realistic values of exchange couplings j as seen in Fig. (JTl). 



CQ 




16.6 16.8 17.0 17.2 17.4 17.6 17.8 18.0 
W 

FIG. 7: (Color online) The quantum phase diagram in a magnetic 
field B (arbitrary units), as function of disorder amplitude W, in 
units of t. We set d = 3, r; = 2 and j — .2. 
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X. FINITE TEMPERATURE PROPERTIES 

We can now proceed to calculate finite temperature proper- 
ties. To this end, we first derive the density of free magnetic 
moments at temperature T. It can be obtained by integrating 
P{Tk) according to 



nFMiT)^nM dTxPiTK) 



(59) 



At low temperatures, T <C T^ the free moments are deter- 
mined by the tail of the distribution. 



A. Insulator 

In the insulating regime for T > A^, P{Tk) is given by 
Eq. ( 25 I, while at smaller temperature, not exceeding the lo- 
cal level spacing, T < A^, the tail of P{Tk) is changing to 
Eq. (36 1. Thus we find the density of free moments in the 



insulating regime 



npAiiT) = nFMiO)+nM 



]{0 r^<*'^forT<A, 



(60) 

where the density of free moments at T = 0, rtpM (0) is given 
by Eq. (|37]l, decaying towards the AMIT. This result is plotted 
as function of disorder amplitude W for various temperatures 
T in Fig. (ISll. Thus, we find that the magnetic susceptibility is 
diverging at low temperature 

,„^ npAiiT) 1 

X{T) = nFAiiO)- + 



riM 



Er 



m 



^-K7(*)'for T < A^ 



j-i 



(61) 



with a Curie tail, whose weight increases as the Fermi energy 
moves deeper into the insulating regime. The specific heat is 
given by. 



Accordingly, we find that the magnetic susceptibility x(r) — 
nFM{T)/T has also in the metallic regime a power law tail, 
which is however cutoff at T < A^, where it converges 
to zero. We get the contribution to the specific heat, using 
(^(-y>) _ y d«FM(T) ^ nFM(T), vanishing at low tempera- 
tures T < A^. 

One may ask if these thermodynamic results are modi- 
fied by inelastic scatterings, introducing a thermal length Lt, 
which decreases with increasing temperature according to 
Lt ^ T^^l^ , where z is the dynamical exponent. If electron- 
electron and electron-phonon scatterings are disregarded, Lt 
is set by the relation T = A^^ = D/L"^^, yielding z — d. 
Therefore, it has been argued that on length scales exceeding 
Lt, the system size L and the localisation length/correlation 
lengths ^ are substituted by Lt, in the scaling theory of the 
AMIT^^. Thus, indeed, for temperatures T > A^, ^ would 
have to be substituted by Lt- However, as we find above, in 
this temperature range the results do not depend on ^ anymore, 
so that the finite Lt does not modify the above results. 
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FIG. 8: (Color online) The ratio of free magnetic moments (magnetic 
moments with Tk < T), npM/nM, as function of disorder strength 
W on the insultating side of the transition for temperatures T/ Ec = 
0, 1 X lO"'^, 5 X 10"'', 1 X 10"^ from the bottom to the top curve. 
We set j = 0.2, d = 2 = 3 and 77/d = 2/3. 
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dT 
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for T > A^ 

(62) 



These results also apply at the AMIT, where npM{0) 
and A^ -^ A. 



B. Metal 

In the metallic regime, P{Tk) is for T < A^ given by Eq. 
S2I1, so that we obtain. 



npM{T) 
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■n 
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exp 
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(63) 



XI. FINITE TEMPERATURE PHASE DIAGRAM: 
KONDO-ANDERSON TRANSITIONS 

The spin relaxation rate is found to depend on tempera- 
ture, Eq. ( [44| , due to the Kondo screening, and the tem- 
perature dependence of the density of free magnetic mo- 
ments, Eq. (60 1. Therefore, the amount by which the spin- 



and time-reversal-symmetry is broken, depends on tempera- 
ture as well. Since the position of the AMIT is determined 
by these symmetries,**"* it shifts as function of temperature 
T: when the spin relaxation rate increases with temperature, 
the AMIT shifts towards larger values of disorder amplitude 
W. Thus, a metal-insulator transition may occur at a finite 
temperature Tc{W, J). In order to investigate the existence 
of such a transition, we will apply the Larkin-Khmel'mtskii 
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conditiorr^^ with the temperature dependent symmetry pa- 
rameter ^^.(r) = S,^ / Df,Ts{T). Since the condition Xs = I 
gives an estimate for the position of the transition, Wc{J, T), 
we find. 



[£,{W,(J,T)f = al 



W^ 



WMT)-WO 



2v 



DeTs{J,T). 



(64) 

Again, as in the previous section, we can apply this criterion 
in two ways: 



A. Approaching the AMIT from the Insulator side. 

Coming from the insulating side of the transition, where the 
localization length ^ is still finite and smaller than the ther- 
mal length Lt, the ratio X = ^"^ / Di,Ts{T) is finite, giving 
a measure of the amount of time reversal symmetry break- 
ing. \/Tg{T) saturates at low temperatures to the spin relax- 
ation rate from free magnetic moments Eq. (J47]i. Thus, at low 
temperatures the transition occurs at W^{J) as given by Eq. 
(49 1. At higher temperatures 1/ts{T) increases with a power 
which depends on the distribution of the Kondo temperature. 
Thus, the spin relaxation rate at finite temperature is given by 
a weighted integral over the distribution function of Tk, as 
1/ts{T) = J^ dTKP{TK)l/Ts{TlTK). The spin relaxation 
rate of a magnetic impurity with a given Kondo temperature 
Tk, 1/ts{T/Tk), follows Eq. ^. Thus, it increases first 
like T^/T^ when T < Tk, until it reaches a maximum and 
then decays logarithmically slowly towards its classical value 
2^^ciassicai gjjj^-g q^j- analysis is limited to low temperatures 
T ^ T'^, we can simplify l/r^ (T) as a sum of spin relaxation 
of free moments of density Eq. ( [60| , at sites whose density of 
states is suppressed as p(r) ~ ^^^"fm^ ^nd the spin relax- 
ation from spins whose Kondo temperature exceeds T, which 
we can approximate by T]^ due to the peaked distribution. 
Thereby we get in good approximation the spin relaxation rate 
in the insulator as 
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(65) 
When Xs — 1 the symmetry breaking is sufficient to shift 
the transition to the larger disorder amplitude Wl{J, T). Ac- 
cordingly, we find the transition temperature. 



T/ = E,ci 



W -W^ 
W? 



(66) 



where W^ is given by Eq. 

{kj/j)~^/^{D^t'^Y^/3, where I/k^ 



(49 1 



17(2 



and c/ 



B. Approaching the AMIT from the Metallic side. 

Coming from the metallic side, the density of free magnetic 
moments is decaying fast at temperatures T < A^ according 




W 



FIG. 9: (Color online) The finite-temperature phase diagram of 
Kondo-Anderson transitions. The solid lines are plots of the crit- 
ical temperatures T^' {W, J), Eq. l[68}, and T^{W, J), Eq. ^ , 
respectively, where the disorder amp itude W is given in units of the 
hopping parameter i and the temperature is in ratios of Ec- We used 
the following parameters: j = 0.2, ao = 4, d = 3, ri/d — 2/3, 
u = 1.57 and a'i/{DeT^s°^) = 0.1. 



to Eq. ( |63) . Thus, the spin relaxation rate is dominated by the 
screened magnetic impurities, yielding. 
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where I/t" is given by Eq. (48 i. Thus, the transition is shifted 
to W^^ (T) and, accordingly, we find a transition temperature. 



T^^'iW) - ^S{S+1)tt^D, 



w -W^ 



wp 






(68) 



This is plotted in Fig. l9]as function of disorder amplitude W . 
Thus, we can conclude that there is a critical semimetal 
region which extends over a finite temperature range, T/ < 
T < T*^. Since we derived the scaling function only at 
small symmetry breaking parameter Xg, the phase diagram 
at larger disorder, where the critical disorder of the unitary 
ensemble (when time-reversal symmetry completely broken) 
is approached, might be modified. This part of the phase di- 
agram is further complicated by the fact that 1/ts reaches a 
maximum and decays logarithmically at temperatures exceed- 
ing T^ . Furthermore, as inelastic scattering and dephasing 
processes will become stronger at higher temperatures, that 
higher temperature part of the phase boundaries is expected 
to become less well defined, mainly indicating a crossover re- 
gion. 



XII. CONCLUSIONS AND DISCUSSION 

We conclude that spin correlations in disordered metals 
do strongly modify the disorder induced Anderson metal- 
insulator transitions. Starting from the numerically well- 
established fact that a system of noninteracting electrons in 
a nonmagnetic disorder potential undergoes a second-order 
transition between a metallic state and an insulator, we studied 
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the mutual influence of Kondo correlations and the Anderson 
localization transition. Since the position of the AMIT de- 
pends very strongly on time-reversal and spin symmetries, we 
find that the critical density and the critical disorder depend on 
the density of magnetic impurities, since these randomize the 
conduction electron spins. However, since the magnetic im- 
purity spin in a metal is screened at low temperatures due to 
the Kondo effect, restoring a renormalized Fermi liquid with 
full time-reversal symmetry,E2 we conclude that coming from 
the metallic side, the transition occurs at the critical disorder 
of a time-reversal invariant system, W^. However, coming 
from the insulating side, it occurs instead at the stronger crit- 
ical disorder of a system where this symmetry is to some de- 
gree broken by free magnetic impurity spins. Since a shift in 
the AMIT results in an exponential change of the finite tem- 
perature resistivity, this has very strong experimental conse- 
quences. Since both the correlation length and the localization 
length are infinite in the critical region in-between the two crit- 
ical points, the zero-temperature conductivity vanishes in this 
region, making the system a semimetal. 

Taking into account the multifractality of the critical wave 
functions, we derived the concentration of unscreened mag- 
netic impurities, the resulting spin scattering rate, and there- 
fore the shift of the AMIT. Information on the multifractal 
distribution of the critical wave functions, which is well es- 
tablished by numerous careful numerical calculations,'^ al- 
lowed us also to characterize this critical region in more de- 
tail. Unscreened magnetic moments occur at sites with low 
intensity, which break the time-reversal symmetry of the con- 
duction electrons. At sites with high intensity local singlets 
are formed, where one conduction electron is captured and 
localised, leaving the symmetry of the other conduction elec- 
trons unchanged. Thus, the Kondo- Anderson transitions share 
some features of a first-order phase transition: In the criti- 
cal semimetal phase there is coexistence of different magnetic 
states of the magnetic impurity spins, which change the states 
of the conduction electrons, accordingly. A magnetic Zee- 
man fields is found to shift the transition to larger disorder W, 
and is dominating the magnetic field dependence for small ex- 
change couplings j < jz- 

An experimental consequence of the critical semimetal 
phase is the divergence of the dielectric susceptibility, x{T = 
0) ~ ^2 _ [VF/(J)-VF]-2'' at the disorder amplitude T^/(J) 
(or at a critical electron density n],{J)), while the zero tem- 
perature limit of the resistivity diverges with the correlation 
length as p{T = 0) - ^c -- (VF - W^)'" at the weaker 
disorder amplitude W^ (or at a larger critical density n^, ac- 
cordingly). This difference, Wl{J) — W^ , increases with the 



concentration of Kondo impurities, as given by Eg. (49 1, un- 
til it converges to its limiting value W^ — VF^.^We expect 
the critical semimetallic phase to be observable in materials 
where both the AMIT and the Kondo effect are present si- 
multaneously at experimentally accessible temperatures, such 
as in amorphous metal-semiconductor alloys^' ^- with dilute 
magnetic impurities,^-' or in doped semiconductors, such as 
Si:P, where thermopower measurements are consistent with 
the presence of Kondo impurities with {T^) ~ IK*^ 

On the insulating side of the transition, physical properties 



such as the finite-temperature resistivity are governed by the 
single-particle gap A^ = l/{pi'^) = E^[{W - W^/W^Y'^, 
which causes activated behavior. 

Aspects of the Kondo-Anderson transitions may be rele- 
vant for the understanding of metal-insulator transitions in 
real materials like Si:P. There have been detailed observa- 
tions of a finite density of magnetic moments in Si:P across 
the transition, from magnetic susceptibility and specific heat 
measurements.^'' These have been previously modeled by the 
phenomenological two-fluid model of Refs. 49 58 which did 
not take into account the Kondo effect. On the other hand, a 
realistic description of Si:P should consider that uncompen- 
sated Si:P contains a random array of half-filled donor levels 
that are randomly coupled. At low-donor concentration, close 
to the MIT, these levels are only weakly hybridized with the 
conduction band. Therefore, Si:P should at least be described 
by a half-filled Hubbard model with random hopping and 
on-site potentials.^i'25|58l'J'JJ^g^ ^jj localized states carry single 
electrons with spin 1/2 and double occupancy is prevented by 
the on-site interaction U . For the Hubbard model without dis- 
order, the ground state at half filling is antiferromagnetic and 
there is an excitation gap of order U . Dynamical mean-field 
theory (DMFT) calculations show that, at least in the limit 
of infinite dimensions, spectral weight is transferred into the 
middle of the gap; it has been argued that these states yield ex- 
tended quasiparticles which may form a Fermi liquid. ^^ These 
quasiparticles interact with the localized magnetic moments 
in the lower Hubbard band, similarly to normal conduction 
electron quasiparticles which interact with magnetic Ander- 
son impurities. In the presence of disorder, it has been shown 
that even when there are no resonant levels, and the donor lev- 
els have merged with the conduction band, the localized states 
in the band tail may carry localized magnetic moments.*"" Cal- 
culations of the disordered Hubbard model, using the statis- 
tical DMFT,'*^' and a variant of the local DMFT methocP 
show also that this model shares many features with the dis- 
ordered Anderson impurity model studied in this article. As 
mentioned above, starting from the Anderson impurity model, 
we neglected the magnetic impurity scalar potential K given 
in Eq. (J9]l. We can justify this now, by referring to nonpertur- 
bative numerical studies of the Kondo effect with a pseudogap 
of power r which find that, below the critical exchange cou- 
pling Jc = rD, a scalar scattering potential is irrelevant and 
the magnetic moment remains unscreened,^^ leaving our con- 
clusion about the density of free moments unchanged. In the 
other case of power-law divergencies, which we find to occur 
in the multifractal states at sites with large intensiy {a < a^), 
a nonperturbative renormalization group analysis shows that 
the on-site potential is irrelevant for all antiferromagnetic ex- 
change coupling amplitudes, ^^ leaving our conclusion on the 
formation of local singlets at such sites unchanged as well. 
At finite density of magnetic moments the indirect exchange 
coupling is in competition with the Kondo screening. In a dis- 
ordered metal these couplings are widely distributed. In the 
vicinity of the AMIT, their distribution function can be calcu- 
lated in a similar manner than the Kondo temperature, and this 
will be presented in the subsequent article Ref. TS" 

A nonphenomenological, analytical approach to the AMIT 
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of interacting, disordered electrons is provided by field theo- 
retical methods that typically employ the replica trick.^"* De- 
coupling the resulting action of interacting Grassmann fields 
with the Hubbard-Stratonovich transformation, a subsequent 
gradient exp ansion of the matrix fields results in the nonlinear 
sigma model, !^ ! ^ . This kind of field theory treats the disor- 
der scattering nonperturbatively and allows for a renormal- 
ization group analysis. It turns out that not only the conduc- 
tance but also a frequency parameter and the spin density are 
renormalized. In particular, it has been found that the disor- 
dered electron system becomes unstable to local spin density 
fluctuations at the AMIT.^'* It is still unclear if these indicate 
the onset of a magnetic phase transition,^* *^ or are rather pre- 
cursors of the formation of local magnetic moments!^ Since 
the nonlinear sigma model is formulated in order to account 
for long wavelength fluctuations such as charge diffusion and 
spin diffusion modes, it cannot, at least in these early formu- 
lations, describe local moment formation. Therefore, it is an 
open challenge, to combine the field-theoretical approach to 
the AMIT with a nonperturbative treatment of the local mo- 
ment formation and the Kondo correlations caused by them. 
Field theoretical formulations of the disordered Kondo prob- 
lem have been developed in the mean field approximation,!^^!^ 
and in the Ginzburg -Landau theory.^" ^' However, in these ap- 
proaches, the AMIT of the conduction electrons in the disor- 
dered potential has not yet been taken into account. We hope 
that the theory of the Kondo-Anderson transitions presented 
in this paper, which was based, to some extent phenomeno- 
logically, on the multifractal distributions and correlations at 
AMITs, can serve as a guide in the quest for a theory, which is 
derived from the Hamiltonian of disordered, interacting elec- 
trons, taking into account nonperturbatively both disorder and 
interaction effects, including Kondo correlations. 



Appendix A: Wave function Correlations in the vicinity of the 
AMIT 



1. One Energy at the mobility edge Em- 

The long-range spectral correlations in a d-dimensional 
system can be quantified by spatially integrating the correla- 
tion function of the eigenfunction probabilities associated to 
two energy levels distant by conm = En — ii'm-'^ When one 
of these energies is at the mobility edge Ej^j, one finds 
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Max(|w„„|,A)' 

rim I J 1 
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Ep > \uj„ 



>E„ 



(Al) 



where < t] < d, with rj = 2{aQ — d). This exponent is 
obtained by the requirement that in the limit of small energy 
differences, |w„m| ^ A, one recovers 



C„™||._|^A->i''(l^n(r)r)^i 



^ d—T2 



(A2) 



where T2 = d2 = d — 2(q;o ^ d) and A = D/U^. For 
\'^n'm\ < Ec, correlations are enhanced in comparison to the 
plane-wave limit, where C„m = 1- Note that for \ujnm\ > Ec 
the correlation function decays below 1. This anticorrelation 
ensures that the intensity is normalised: A dip in the intensity 
at one energy implies an enhancement of intensity at another 
energy in the band. This anticorrelation is expected to occur 
up to some finite energy Ep, beyond which the correlation 
function increases to the uncorrected value 1 in a nonuniver- 
sal way. 



2. Both energies at a distance from the mobility edge E^ 
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The spectral correlations in a d-dimensional system also ex- 
ist away from the transition whenever either the correlation 
length ^ (on the metallic side) or the localization length ^^ (on 
the insulator side) are finite. In the correlation function of Eq. 
( Al I, the energy difference \u!nm\ is for e„ < e,„ then substi- 
tuted by Max[|w„m|, A^^], where A^^ — Ec {^ji/a.c)^'^, and 
Uc is defined by Ec ~ l/(/5aJ?), with p denoting the average 
density of states. We get therefore 

Cnrn = L'' J d'r {\Mr)\^\^^m{r)\^) 



Eg 



Max[|w„„,|,Aj^] 
\ 2 

Eg 



, < \uJnm\ < Ec, 



(A3) 



IWnml > Ec 



where rj = 2{aQ — d). For conm < Ec, correlations are 
still enhanced in comparison to the plane-wave limit, where 
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3. Joint Distribution Function. 

The joint distribution function of two eigenfunction inten- 
sities at the same position should be log-normal since the dis- 
tribution function of a single eigenfunction intensity is log- 
normal as given by Eq. (|2]i. Therefore, we recently conjec- 
tured it recently to be of a log-normal form, when one of the 
energies is at the mobihty edge.^' In general, when, one or 
both energies are away from the mobility edge, the finite cor- 
relation length/ localisation length ^ needs to be taken into ac- 
count when one or both energies are on the metallic/localised 
side of the transition, respectively. We can then conjecture the 
joint distribution function for a; and ak, where, in the metal- 
lic regime, we define a — ~ "*^ ', l^^'^'' ' which has the 
'=' ' 1114 ' 

distribution P{a) ^ ^-(a-ao)V(2»j) Accordingly, in the lo- 
calized regime, we get the same distribution function, defining 
there a ~ — ln[|V'(a;)p]/ln^. Then, for |e; — ek\ < Ec, we 
conjecture the joint distribution function to be of the form. 

It/ \ J1 \tl \ Jl — (7,, (°l-°0)(°fc-°0) 

r-yai.akj — C,i ?j, ^n^ , 

(A4) 
where ^/ is the correlation/localization length of a state at en- 
ergy El, Kik - Max [\Ei - Sfc|,Min(A4,, AjJ] /E,, and 
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We consider first the case, where, one of the energies is fixed 
to the mobility edge, En^ — Em- The power law of the corre- 
lation of higher moments of the intensity, can then be obtained 
by taking the limit of small energy differences, a;„._„ — > A, 
which yields. 






~ l(?-i)'*~^' -L''". (AlO) 



Tq is found to terminate for q > qc = ao/r]. Thus, for d ~ 3 
and ao ~ A, qc ~ 2, so that all higher order moments have the 
same Tq ~ T2 = d — ri,yq > qc- Thus we find 



r]q = {q-2)d + ri. 



(All) 



Furthermore, we notice that the pair of intensities whose en- 
ergy are closest to each other are correlated with a power rj/d. 
Thus, we can conclude that the 3-rd order correlation function 
is to leading order given by 
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This result is valid for < |w„; „ | < E^i,j. If all ener- 
gies in C^ are placed away from the mobility edge, En^ ^ 
Em, then the energy difference |Min(aj„..„ ) |, is replaced by 
A„. = Ec/£,n. if it is smaller than that energy scale, respec- 
tively. 



where we introduced the length scale ^ through 
ln|= Vln^lnCfc- 



(A6) 



Averaging the local intensity of the state with energy Ei 
with the conditional probability P{a,ai)/P{a), which can 
be rewritten as 



P. 



Ml) =Ci 



dlnd 



= 0)1" 



we then get Eq. (28 1 
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Note that when the energy Ek is at the mobility edge 
Ek = Em then f j. — > L, in Eq. ( |A4[ ), and we recover the con- 
ditional intensity at energy Ei, Eq. (jSl), in a more rigourous 
way (the finite correlation length of that state, ^i, did not ap- 
pear explicitly in our previous conjecture Eq. (7) of Ref. |2T1 
but is now taken into account in Eq. (|A4|i). 



4. Higher Moment Correlation Functions 

The correlation of higher moments of the intensities is de- 
fined by 



^u = i'^'dii^-wi' 



(A9) 



Appendix B: F [a, Tk] 

Let us first consider the case, when the Fermi energy is 
at the mobility edge. We can then evaluate F[a,TK] in 
the following way. We transform the summation over states 
I to the integration over energy e. Then, we transform to 
t = —\n{e/Ec)- Next, we can substitute in good approxi- 
mation tanh X K, X for < x < 1, and tanh a; « 1 for a; > 1. 
Thereby, we get 
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X exp[— 
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The integrals can now be performed, yielding Error functions. 
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ErflVAa + B/{2VA))]), for each term in Eq. (|Bl). Since 
A = {ao — d)/{d'^ In L) the arguments of all Error functions 
diverge for L —>■ 00, and we can use the asymptotic expansion 
Erf[z] — >• signz x 1 — e~^ /{^/nz). Thereby, we find for 
a > ao — d. 
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where da ^ 
then for L 



d/{ri + 2a). Solving 1 = F[a,TK] we obatin 
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